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We have analyzed and designed an organic high 
speed digital optical switch (DOS) based on trans- 
verse electro-optic effect. In analysis section, we pro- 
posed a quantum photonic model (QPM) to explain 
linear electro-optic (EO) effect. This model inter- 
pret this effect by photon-electron interaction in at- 
tosecond regime. We simulate applied electric field 
on molecule and crystal by Monte-Carlo method in 
time domain. We show how a waveguide response 
to an optical signal with different wavelengths when 
a transverse electric field applied to the waveguide. 
In design section, we configure conceptually a 2 x 2 
EO switch with full adiabatic coupler. In this DOS, 
we use a rib waveguides that its core has been con- 
structed from NPP crystal. This switch is smaller at 
least to one-half of similar DOS in dimensions. 
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High-speed and high-performance optical switches 
are vital for optical communication networks and op- 
tical signal-processing systems. Among these, op- 
tical switches based on the electro-optic(EO) effect 
are promising for applications that require high-speed 
low- loss switching such as: optical burst switch- 
ing, optical packet switching or EO beam deflec- 
tor switching [l]-[16]. MachZehnder interferometer 
(MZI) switches and directional coupler(DC) switches 
have been under intensive investigation for a long 
time, but they both have very limited fabrication tol- 
erance. Furthermore, bias control networks and pre- 
cise switching voltages are needed in MZI switches 
and DC switches to achieve high performance p]lfT4]. 
Electro-optic polymers are particularly interesting for 
new device design and high-speed operation [lj-[7j. 
Organic optical materials like MNA, NPP, MAP have 
a high figure of merit in optical properties in compar- 
ison with inorganic optical materials such as BBO, 
LiNbOa, [17J. 2-methyl-4-nitroaniline (MNA) and 
N-(4-nitrophenyl)-L-prolinol (NPP) have the highest 
figure of merit between organic nonlinear optical ma- 
terials [18]. Thus they are used for electro-optic and 
nonlinear optic applications [T7| - [3T] . Because of the 
large electro-optic coefficient of organic material, a 
certain amount of refractive index (RI) change can 
be realized with lower driving voltage than in other 
EO materials(1.51///im for NPP,[19j). Since there is 
better optical/RF velocity matching inside the EO 
organic, high-speed switching can be achieved eas- 
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ily. Digital optical switches (DOSs) are based on 
adiabatic propagation or modal evolution effect [8j. 
These switches are wavelength-polarization insensi- 
tivity and electronics for control of these devices are 
simple. These switches can be applicable in optical 
network of phase array radar systems. Their struc- 
ture of waveguides can be chosen the buried-type 
waveguide or rib- type (ridge- type) waveguide. In this 
paper we analyze and design an organic DOS that 
is based on EO phenomenon. Our favorite organic 
molecule is NPP. In analysis, we use QPM, to ex- 
plain linear optical phenomena in molecular scales. 
We show that the phase retardation of input light 
with different wavelengths is distinctive, when it trav- 
els through the waveguide. In design we design a 2 x 2 
DOS with adiabatic couplers. The main advantages 
of this switch are wide optical bandwidth, low power 
consumption and high speed data transmission. 

2 Analysis 

Some approaches are used to analyze and design op- 
tical switches in micrometer and nanometer sizes by 
quantum mechanics (QM) [32]. In this paper an EO 
switch has been analyzed using QPM. This model 
gives us a constitutive vision about phenomena and 
real materials [33] -[36]. This approach is based on 
four elements: 1-Quasi-classic principle for justifica- 
tion of optical phenomenon in molecular scales, 2- 
Knowledge of crystal network and its space shape, 3- 
Short range intramolecular and intermolecular forces, 
4- Monte-Carlo time domain simulation. In this 
model we suppose a laser beam is a flow of pho- 
tons while passing through single crystal film, inter- 
acts with delocalization 7r-electron system of organic 
molecule and delays the photon in every layer [36] . 
By precise calculation of these retardation in every 
layer, we obtain RI in specific wavelength and ex- 
plain EO effect too. After one, we apply an external 
electric field to crystal to change its RI a little to 
switch specific wavelength to identified channel. 

2.1 Crystal and Molecular Structure of our 
Favorite Organic Compound 

NPP (Ciii^i4A^203) crystallizes in the solid 
state in an acentric monoclinic (with space 
group P2i) structure and their parameters are: 
a=5.261yl°,b=14.908A°,c=7.185A°, /3=105.18° and 
in the wavelength range of 0.5 to 2 jam is transparent. 
The most interesting property of NPP crystal is 




Fi gure ll The crystal packing of NPP. This molecule has a 
monoclinic crystal structure with space group P2i. There are 
two molecules in every cell. 
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Figure 3: The NPP molecule in dielectric frame. The bond 
lengths (in angstrom) and angles of NPP molecule have been 
shown too. As we see the angle between charge transfer axis 
and Y axis is 58.6°. The mean plane of molecule is in X-Y 
plane approximately (11° accurately). 

the proximity of the mean plane of molecule with 
the crystallographic plane (101); the angle between 
both of these planes being 11°. Nitro group of 
one molecule in downward connects to Prolinol 
group in upper by hydrogen bonding. The angle 
between b orientation of crystal and N(l)-N(2) axis 
(charge-transfer axis) is equal to 58.6° [17]. FigH] 
shows the molecule and crystal packing of the NPP. 
For accurate and valid simulation, these properties 
and angles have to be exerted. For benzene molecule, 
benzene ring is a circle (see figEJ; In NPP molecule, 
nitro group acts as an acceptor and the other main 
groups on the other side of benzene ring acts as a 
week donor (see Fig(T]). Fig l3] demonstrates the bond 
lengths and angles of the NPP molecule. As we 
see in this figure the bond lengths in benzene ring 
are not same. In our simulation, for similarity we 



2 




(a) 




(c) 




Superposition of (a), (b) and (c) 

Figure 2: Electron cloud for Benzene molecule that obtained from Huckle theory [37] and its approximation. 



consider an ellipse correspond to circle for electron 
cloud. We obtained £=0.26 for the ellipse of NPP 
from simulation. 

2.2 Our Model for Electro-Optic Effect 

Assume that single crystal film lies in x-y plane and 
light propagation is in the z direction; therefore in 
the presence of an electric field the equation of index 
ellipsoid by assuming crystal symmetry will become: 

[(^).+ri2^^]x2 + [(-^2),+r22^^]^' + 2r6i^.xy = 1 

(1) 

that Ex and Ey are transverse electric field com- 
ponents, [19j. With appropriate rotational trans- 
formation, these relations can been simplified. In 
NPP, ri2 and r22 are large coefficients. There- 
fore Tlx and Tiy simplified to Ux — (\l2)n\r\2Ey and 
Tiy — {l/2)nyr22Ey respectively. The phase retarda- 
tion F, with an applied electric field in a typical linear 
transverse electro-optic modulator will be obtained as 
follows, [38i: 



F = 



1 



{nlr22-nln2)Ey] (2) 



The total phase difference between two perpendicular 
polarization of light (in our example Ex and Ey)^is 

A = Ao + (5 (3) 

where Aq is due to linear birefringence and S is due 
to linear electro-optic effect. In this case S is much 
smaller than Aq. 

In sub-micron space scales and sub-femtoseconds 
time scales, the optical constants loses its stabiliza- 
tion and classical equations for linear and nonlinear 



optical phenomena are not useful [36j. 
Now, we suggest a microscopic model for linear 
electro-optic phenomenon. That would be assumed 
only linear optic phenomenon exist and nonlinear op- 
tic phenomenon do not exist, approximately. Because 
laser watt is not much, one photon interacts with any 
NPP molecules. In each interaction between photon 
and electron in every layer of crystal, we suppose a 
delay time equal to {ith layer of the crystal). Total 
delay time for m layers in crystal region is equal to: 



Consequently required time for photon transmission 
in L length of crystal is equal to r, achieved from 
relation: 



r = 



nL 

Co 



L 

Co 



E' 



(4) 



where Cq is velocity of light in vacuum. By using 
this relation, we can relate macroscopic quantity n 
to microscopic quantity r^. In biaxial crystals, r and 
consequently n depends on polarization direction of 
incident light. Because dipole-field interaction con- 
clusion is different for any direction of molecule. If 
Tx would be a microscopic delay for interaction of 
x-polarization field with dipole (or charge transfer 
action) and Ty would be a microscopic delay for in- 
teraction of ^-polarization field with dipole then the 
final phase difference between these two fields (named 
phase retardation) will be: 

A(l) = uj'{rx-ry) =uj- ^(r^, - r^J (5) 
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Of course this relation give Aq of (6). We justify S 
from our model in later subsections. 

2.3 Photon-Electron Interaction in Attosec- 
ond Regime 

It can be shown that for absorbtion or emission of 
photons the material has to perform a transition be- 
tween two eigenstates Em and Ep of the material and 
thus the photon energy Ephoton has to fulfill the res- 
onance condition: 

-^photon hUpiioton ~ \^p ^m\ (6) 

But for our linear phenomenon, the photon energy 
is about 2eV (in A = 630nm). If electron would 
be in HOMO (Highest Occupied Molecular Orbital), 
this electron do not go to LUMO (Lowest Unoc- 
cupied Molecular Orbital) or exited state by in- 
teraction. This phenomena is named nonresonant 
phenomenon, [iHl [40l [41] (nonresonant phenomena 
is not exclusive for nonlinear optical phenomena). 
Therefor electron after interaction, go to quasi states 
that their life times is very short, then this electron 
go back to primary state after very short time. The 
nonresonant lifetime is determined by the uncertainty 
principle and the energy mismatch between photon 
energy in second time and the input photon energy. 
We can assume that the characteristic response time 
of this process is the time required for the electron 
cloud to become distorted in response to an applied 
optical field. This response time can be estimated as 
the orbital period of the electron in its motion around 
the nucleus which is about r 10~^^s or 100as,[41j. 
We can estimate this characteristic response time ac- 
cording to (4) if, n = 3, L = 3/i. Consequently Sr is 
equal to 10~^^sec. Because in b direction of crystal 
in Sjnm length, approximately 4024 molecules exist, 
therefore the average quantity of r: 



is in order of 10~^^5 or 1 as. The perturbation in this 
very short time can assume semiclassically. In lin- 
ear phenomenon in this short time, just one photon 
interacts with one molecule. Because NPP molecule 
has delocalization electrons, (or 7r-electron system), in 
benzene ring, that photon interacts with this elec- 
tron type, 02] and it is annihilated [38j. We call this 
photon, a successful photon, (that does not produce 
phonon). 




Figure 4: Assumed 7r-electron orbit of NPP molecule that 
obtained from Fig|2]for Benzene molecule (approximately). 

To obtain 7r-electron wavefunction for benzene 
molecule the Schrodinger equation may be solved. 
Since this is very complicated process, it cannot be 
done exactly, an approximated procedure known as 
Hiickle method must be employed. In this method, 
by using Hiickle Molecular-Orbital (HMO) calcula- 
tion, a wave function is formulated that is a linear 
combination of the atomic orbitals (LCAO) that have 
overlapped [37] (see Fig. 3): 

^ = Y,C^^i (7) 

i 

where the refers to atomic orbitals of carbon atoms 
in the ring and the summation is over the six C atoms. 
The is the probability of the 7r-electron at ith 

atom. Thus: 

|Ci|2 + |C2|2 + \Cs\' + |C4|2 + + ICel^ = 1 

In the case of Benzene molecule: 

as followed from the symmetry of the ring [44l |45] . 
But NPP molecule isn't such as Benzene molecule. 
NPP is polar molecule. Nitro {NO2) is more power- 
ful electronegative compound than prolinol and pulls 
TT-electron system; consequently, the probability of 
finding 7r-electron system at various carbon atoms of 
main ring isn't the same and the probability of find- 
ing TT-electrons near the Nitro group is greater than 
near the prolinol group. Therefore there is no symme- 
try for NPP and electron cloud is spindly or oblong, 
(similar to dom-bell) (FigHj). We estimate this form 
of electron cloud by an ellipse that our calculations 
would be uncomplicated. We assume effective posi- 
tive charge that is located in one of focal points of 
ellipse. The quantity of this effective positive charge 
is determined by semiclassical arguments. For at- 
taining probability of electron presence on an orbit 
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Figure 5: 7r-electron system approximation for NPP 
molecule, ^ = is in the positive direction of dipole (Proli- 
nol side) and ^ = tt is in the negative direction of dipole (Nitro 
side). 

(FigE]), we say, T time is required by radial vector 
to sweep total tt.u.v interior area of ellipse (u and v 
are semimajor and semiminor axis of ellipse respec- 
tively), in t times, this radial vector sweeps: 

t 

TT.U.V.— 

T 

area of ellipse, (see Figl5j). If t is the time, that elec- 
tron sweeps 9 radian of orbit then t is obtained from 
this relation [33j: 

t = — {2arctan(W- tan(-)) 77;^^} 

27r^ Wl + e 2 l + £.cos(l9) ^ 

(8) 

Where e is ellipse eccentricity. By using this rela- 
tion, we attain the required time (t) for electron to 
traverse from 6 to 6 -\- d6 and it is divided by to- 
tal time T. By this approach, we can determine the 
PDF (Probability Density Function) approximately. 
The PDF in apogee (near the Nitro group), is max- 
imum and in perigee (near the Prolinol or Methane) 
is minimum. Therefore PDF is correlated to from 
(8) and seen in Fig|6l The r quantity is correlated to 
and consequently, r quantity is correlated to pres- 
ence probability of 7r-electron system. 
The angle between Y vector and charge transfer ac- 
tion (N1-N2) is 58.6° and X and Z axis is perpendic- 
ular to Y (Fig 13]). We consider propagation along Z 
direction. We spot a photon interacts with 7r-electron 
of NPP in first layer. After interaction, this photon 
gives its energy to electron and is annihilated. Elec- 
tron absorbs energy and digresses in direction of pho- 
ton momentum. Electron with photon energy, may 
not be unbounded and after arriving to apogee of di- 
gression, it returns back to ground state, because the 
photon energy is equal to hu = 1.96ev (h is Planck's 



oh 

1 2 3 4 5 6 

Teta in Radian 

Figure 6: PDF of presence probability of 7r-electron in dif- 
ferent energy states when it rotates around the assumed orbit 
in FiglSl as we see the PDF is correlated to 0. 

constant and u is frequency of laser beam) whereas 
energy for excitation is greater than 3ev. When elec- 
tron returns to ground state one photon is produced. 
The time coming up and down is r delay time. This 
photon after freedom goes to second layer in direc- 
tion of annihilated photon (nonce, we assume the po- 
larization doesn't change), in second layer this pho- 
ton interacts with another delocalization 7r-electron 
certainly, because the effective interaction range of 
photon is approximately equal to its wavelength and 
is very greater than the distance between molecules. 
This molecule is nearest to photon effective central. 
This action is repeated for each layer. The location 
of photon-electron interaction is significant in every 
molecule and it is effective on r quantity directly. We 
assume that interacting photon has circle polarization 
and electron subject to virtual positive charge center. 
The phase retardation between and Ey (2) can be 
obtain from this relation. [33j: 

^ , ij.V2hu.m , . M 9 

= KZe^ E[^o«(^^) - sm(^i)]-^-| (9) 

i=l 

that 

_ {l-e^).u 
~ l + e.cos((9i) 

where e is the elliptical eccentricity and u is the 
semimajor axis of the ellipse. By applying an ex- 
ternal transverse electric field to organic crystal (in 
the range of several volts per micron) the shape of 
TT-electron system will be deformed slightly and we 
would be expect some noticeable variations in micro- 
scopic delay parameters (rx^Ty) and phase retarda- 
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Figure 7: Phase retardation between Ex and Ey of optical signal with different wavelength by applied external electric 
field. As we see; by applied electric field, the virtual ellipse eccentricity change (electron in electrical field). With eccentricity 
variation, the microscopic delay (rx^Ty) change and consequently phase retardation is modified. 



tion; (see FigCj). Consequently, by step-like change 
of input voltage, optical signal is switched between 
output port of DOS. 

2.4 Variation Analysis of Phase Retardation 
versus Applied Electric Field 

We simulate phase retardation of 3/im-length NPP 
crystal by Monte-Carlo method, then we generate 
random number using MAT LAB program. This pro- 
gram produces PDF quantities was explained in be- 
fore subsection and relates each of them to every 
molecule. These values are indexing 7r-electron po- 
sitions in each layer, by assumption a reference point 
(see FigEj). Additionally we have used a MAT LAB 
program for Monte-Carlo simulation. The inputs of 
this program are: 

1. The wavelength of incident optical beam in 
which we want to design DOSs; 

2. /i, m, A: = j^^Cq that are Planck's constant, 
electron rest mass, elementary charge. Coulomb con- 
stant and speed of light respectively. 

3. Unit cell parameters of NPP crystal: a, b, c, (3 
and its other parameters that have given in subsec- 
tion (A). 

4. L: crystal thickness that in our simulation it is 
3/im. 

And the outputs of MATLAB program are: phase 



retardation in each wavelength. 

System calibration is done semiclassically by ex- 
perimental refractive index data.t In this method that 
we obtain three refractive indexes with x-polarization 
in threea with e (eccentricity), u (semimajor axis of 
ellipse) and Z (equivalent positive charge) in a way 
that refractive indexes in three wavelengths are very 
close to experimental data. Then we would see that 
refractive index in other wavelengths and other po- 
larization with same e, u and Z will be achieved. 
Of course these values, e, u and Z would be close 
to experimental structure of crystal, for example u 
would be greater than and smaller than minimum 
and maximum sizes of six lengths of benzene hexago- 
nal respectively, or e would be small but greater than 
zero. In other hand these values must be logical. 
From this method in our simulation we have obtained 
e = 0.26, Z = 3.9, u = 1.4A^ that is very close to ex- 
perimental and structural data. 

Xu and co-workers [T9^, ^ , have done some electro- 
optic experiment about single crystal film of NPP. 
They have obtained |n^ri2 — n^r22| = 340pm/V and 
ri2 = 65pm /V in an optical beam with 1064nm wave- 
length. They have studied phase retardation between 
Ex and Ey of optical beam as a function of angle 
between electric field and charge transfer action of 
NPP. They have concluded that the maximum phase 
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Figure 8: The typical waveguide that is used for our switch. 
The core of this waveguide is NPP and the cladding of it, is the 
poly(tetrafluoroethylene) (PTFE)[27 . We assume input light 
that has xy-polarization and travel in zy-plane. 



retardation was observed for the field oriented along 
the charge-transfer axis which was parallel to the film 
surface. The electro-optic effect or phase retardation 
was negligible when the electric field was applied per- 
pendicular to the charge-transfer axis. This concept 
could be justified by our model in previous subsec- 
tion. When the angle between charge-transfer action 
and external electric field change, the ellipse eccen- 
tricity modify and consequently, the phase retarda- 
tion alter. Obviously, from FigHlwhen the external 
electric is parallel to charge-transfer axis, the ellipse 
drag more and the ellipse convert to a line. Therefor, 
ellipse eccentricity arise and from Fig 111 the phase re- 
tardation growth. In the other hand, from Fig. 5 when 
the angle between external electric field and charge- 
transfer axis change, the ellipse is gathered and con- 
vert to circle and the eccentricity decrease to zero. 
Thus from Fig|7l the phase retardation is lowered. 

3 Design 

In this section, we design an organic 2x2 DOS 
that is based on EO phenomenon with adiabatic cou- 
pler. The adiabatic coupler has some advantages such 
as: a more evenly split power can be achieved in 
the switch and comparably longer electrodes can be 
used for the same device length, which leads to ei- 
ther lower driving voltage or more compact device 
design, [3J. The basic theory of this switch was ex- 
plained in previous section. In this switch we use a 
waveguide that has been shown in FiglHl The core of 
this waveguide is NPP and the cladding of it, is the 
poly(tetrafluoroethylene) (PTFE),[27J. We assume 




Green light Red light 



Fi gure 9 1 The layout of our 2x2 DOS with a full adiabatic 
coupler. Driving electrodes are on the top of the waveguides in 
the interaction region. By applying voltage on the electrodes, 
we can change the refractive index in one or both waveguides. 



input light that has xy-polarization and travel in zy- 
plane according to FiglHl 

In Figini the layout of our 2x2 DOS with a full adi- 
abatic coupler has been shown. The input and output 
waveguides do not have been shown. Driving elec- 
trodes are on the top of the waveguides in the inter- 
action region. By applying voltage on the electrodes, 
we can change the refractive index in one or both 
waveguides. If the refractive index difference between 
the two arms in one side is large enough, the power 
will output from the waveguide with higher refractive 
index. For these EO devices, the material is poled in 
the vertical direction (see Fig fTQ|) . Consequently, the 
driving electric field should be in the vertical direc- 
tion, which is from the top driving electrode to the 
bottom electrode, to have a larger EO effect. To en- 
sure the electric field has a larger component in the 
vertical direction than the horizontal, the spacing be- 
tween the electrodes should be no less than the thick- 
ness of the device. The spacings in our devices are 
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Figure 10: Cross section of the interaction region of our 
proposed DOS device. The driving electric field should be in 
the vertical direction, which is from the top driving electrode 
to the bottom electrode, to have a larger EO effect. 

about 10/im. Because the NPP material has large 
linear electro-optic effect, (ri2 = 65pm/ V, that is 
two times of linear electro-optic effect of nonorganic 
material) by smaller change of applied voltage, the 
optical power can switch to the neighbor waveguide, 
or we can shorten the electrode length [4J. Therefore 
we would have expected that our switch is smaller 
than a switch with nonorganic material. The full 
angle between the two neighbor waveguides, AO is 
about Imrad. From Fig|7j we see that the phase re- 
tardation change for different wavelengths, with same 
applied electric field. Therefore we can switch two 
optical beam with different wavelengths to two out- 
put waveguides. By suitable choice of step-like volt- 
age(using an analogue to digital converter(ADC) and 
choosing two of four electrodes), we can switch spe- 
cific beam with optional wavelength to specific output 
channel. 

The cross section of the EO DOS in the interac- 
tion region is shown in FigfTOl As we see, the rib 
waveguide is used for this switch. 

4 Conclusion 

We analyzed and designed a novel 2 x 2 DOS that is 
based on transverse EO effect. In analysis we justi- 
fied linear EO phenomenon by QPM. This suggested 
physical model could be a powerful tool for analyzing 
and explaining processes that happen in waveguides 
with microscopic and nanoscopic sizes. We showed 
how the phase retardation between different argu- 
ments of an optical field with distinctive wavelengths 
can take place. In design, we configured a DOS with 
full adiabatic coupler. We used a ridge waveguide 
that has been constructed from NPP crystal as core. 
By using an ADC and applying four proper voltages 



on four electrodes, one could switch an input opti- 
cal beam with specific wavelength to certain output 
channel. 
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